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In this talk I gave a brief summary of leading order, next-to-leading order and shower 
calculations. I discussed the main ideas and approximations of the shower algorithms 
and the related matching schemes. I tried to focus on QCD issues and open questions 
instead of making a inventory of the existing programs. 

1 Fix order calculations 
1.1 Born level calculations 

The simplest calculation what one can do is the Born level fix order calculation. This 
calculation involves the phase space integral of the tree level matrix element square and the 
jet measurement function. The structure of the cross section is 



where dr^"^\{p}m) is the phase space integral measure, Ai{{p}m) represents the m-paraton 
tree level matrix element and Fj({p}m) is the jet measurement function that defines the 
physical observable. 

This calculation is relatively simple. The integral free from the infrared and ultraviolet 
singularities. The matrix element is basically a complicated expression but it can be gener- 
ated in a automated way. Several implementations can be found in the literature, Alpgen, 
Grace, Helac, Madgraph and Sherpa[2]. 

We can say that the tree level cross sections can predict the shape of the cross sections but 
in general they have several defects: i) Since it is the leading order term in the strong coupling 
expansion the result strongly depends on the unphysical renormalization and factorization 
scheme, ii) The exclusive physical quantities suffer on large logarithms. In the phase space 
regions where these logarithms are dominant the predictions are unreliable, iii) In the 
Born level calculations every jet is represented by a single parton, thus wc don't have any 
information about the jet inner structure, iv) On the other hand in a real measurement, in 
the detector we can see hadrons and every jet consists many of them. We are not able to 
consider hadroniziation effects in the Born level calculations. 

1.2 Next-to-leading order calculations 

We can increase the precision of our theory (QCD) prediction by calculating the next term 
in the perturbative expansion, the next-to-leading order correction (NLO). However this 
is just one order higher to the Born cross section but the complexity of the calculations 
increases enormously. We have to face to algebraic and analytic complexity. 
The naive structure of the NLO calculation is 
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Here a , a and a correspond to the Born, real and virtual contributions, respectively. 
This expression is well defined only in d = 4 — 2e dimension because both the real and virtual 
terms are singular separately in d = 4 dimension, but their sum is finite. Thus we cannot 
calculate them separately, first we have to regularize this integral. In the real part the 
singularities comes from the phase space integral from the regions where a gluon becomes 
soft of two partons become coUinear and the integral over these degenerated phase space 
regions leads to contributions those are proportional to 1/e and . The infrared singularity 
structure of the virtual contributions is exactly the same but with opposite sign, thus they 
cancel each other. To achieve this cancellation we have to reorganize our calculation in such 
a way that can be carried out in d = 4 dimension 

^TNLO= / da^+ / [d^" -da'']^^^+ f / da""] . (3) 

JN JN+1 JN Jl 

Here we subtracted the approximated version of real contribution and added it back in 
different form. In the second term da^ cancels the singularities of da^ and it is safe to 
perform the integral in d = 4 dimension while in the third term the explicit singularities 
of da^ are cancelled by J^^ da^ , where we performed the integral over the unresolved phase 
space analytically. It is important that the approximated real contribution has universal 
structure. This term is based on the soft and coUinear factorization property of the QCD 
matrix elements. A general subtraction scheme was defined by Catani and Seymour and 
the extenion of this method for massive fermions is also available [4 . 

The NLO calculation can be carried out but it hasn't been automated like the Born level 
calculations. The most complicated processes what we can calculate are 2 — > 3 type [5]. 
To go beyond this limit we have to find an efficient way to compute the virtual correction. 
Recently we have had some very promising developement on this area [B] . 

With the NLO corrections we can significantly reduce the dependence on the renormal- 
ization and factorization scales but in some cases it is not enough and the NNLO is also 
required. In these calculations one of the jet is represented by two partons. This can give 
some minimal information about the inner jet structure but is still very poor. The ex- 
clusive quantities are still suffers on large logarithms and we are still not able to consider 
hadronization effects. 

1.3 Next-to-next-to-leading order calculation 

For some processes and/or jet observables it is important to know the cross sections at next- 
to-next-toleading order level. In this cases the NLO if-factor usually large even larger than 2 
which means that the NLO correction doesn't reduce the scale dependences. Recently some 
simple but important processes have been calculated using sector decomposition method (7] 
and there are some ongoing developments on defining a general scheme for NNLO calcula- 
tions [Hj. 

2 Leading order parton shower 

The fix order calculations are systematically defined order by order and usually give good 
description well the data over the phase space where the large px event are the dominant. 
In any order we still have to deal with the presence of the large logarithms and we cannot 
consider hadronization effect. 
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Real time picture 



Shower time picture 



Figure 1: The left-hand picture depicts quark-quark scattering, with time proceeding from 
left to right. The hardest interactions are those toward the center of the picture. These 
are treated first in a parton shower Monte Carlo program. Thus in Monte Carlo time t, we 
start with the hard process and work toward softer interactions (with some of the particles 
moving backwards in real physical time), as depicted in the right hand picture. The rounded 
rectangles represent intervals of Monte Carlo time in which nothing happens. 



There is an other way to calculate crass section in the perturbative framework, the parton 
shower calculations. Consider the parton shower picture of hadron-hadron scattering in 
which there is some sort of hard event, say jet production. The parton shower description 
starts form hard scattering and proceeds forward to the softer scattering. In the final state 
the shower proceed forward in real time but for initial state parton the sowers proceeds 
backward in real time. This is depicted in Figure [T| 

2.1 Shower evolution 

The parton shower evolution can be represented by an evolution equation and it is the 
solution of the following integral equation 

l^{tf,t2)\M2)^m!,t2)\M2)+ dhUiti,t3)nih)Uit3,t2)\M2) ■ (4) 

The shower evolution starts form the hard scattering and it is represented by the function 
|A^2) that is a probability of a given partonic state in shower time ^2- Then U{tf,t2) is 
the probability function of having a particular partonic state in a later evolution time t{. 
The evolution operator is sum of two terms. The first term in Eq. Q represents parton 
evolution without splitting. The non splitting operator JV(t{, 12) that inserts Sudakov factors 
giving the probability that nothing happens between time t2 and t^. The Sudakov is the 
exponentiated inclusive (summed over spin and color and integrated over the momenta of 
unresolved partons) splitting kernel. The second term in Eq. ^ represent the splitting. The 
partonic state is evolved without splitting to an intermediate time ^3 and splitting happens 
given by the splitting operator Ti^t^) and the system is evolved with possible splitting from 
t'i to tf. The splitting operator is based on the universal soft and coUinear factorization 
property of the QCD matrix element. This evolution equation is depicted in Figure |2j 
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Figure 2: Evolution equation. The shower (orange rectangle) starts from the hard matrix 
elements (green rounded rectangle) and the partons are evolved to the final scale without 
splitting (yellow rounded rectengle) or with splitting at an intermediate time (red circle) and 
evolved to the finial scale with possible splittings. 

2.2 Splitting operator 

The splitting operator of the leading order (LO) shower is derived from the factorization 
property of the QCD matrix elements in the soft and coUinear limits. This factorization 
property is universal. Let us start with coUinear factorization. 

When two parton become coUinear the m + 1 tree level matrix element (and the matrix 
element square) can be written as a convolution of m parton hard matrix element and a 
universal singular factor in the spin space. This factorization is depicted in Figure |3] at 
matrix element level. 




Figure 3: CoUinear limit. The universal singular part is represented by the red circle and 
it is based on the 1 — > 2 matrix element at LO level and the yellow ovals represent the 
corresponding m and m + 1 parton matrix elements. 

The coUinear limit has some nice features. In the squared matrix element the singularity 
doesn't make color connection, it is completely factorized out. There are some spin correla- 
tion in the gluon splitting but this spin correlation is rather trivial. On the other hand one 
can always use spin averaged splitting functions but this is an additional approximation. 
Considering only the coUinear emission our first "candidate" for the splitting kerneQ would 
be 

m 

,Pni+l)[Vlm+l(Pl,Pni+l)Y ■ (5) 

1=1 

^Tor the shake of the simplicity I give a formal definition for the spin and color averaged splitting function 
which actually appears in the Sudakov exponent. 
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Figure 4: Soft factorization. The soft gluon is represented by the red wiggle hne and the 
tree level matrix elements by the yellow ovals. The soft gluon can be emitted from any hard 
parton thus it makes color connections all possible way. 



Here T; is the color charge operator of the mother parton and Vim+i is the vertex function. 
Every parton can split and after the splitting the daughter partons are labeled by I and 
TO + 1. The momenta of the daughter partons are pi and Pm+i- With coUinear splitting 
operator one can sum up the leading (double) logarithmic contributions properly. 

The matrix element has universal factorization property in the soft limit, when the energy 
of a final state gluon becomes zero. We cannot neglect soft gluon contributions since they 
produce next-to-leading logarithms when we integrate the matrix elements over the phase 
space. The structure of the soft gluon radiation is depicted in Figure |4] In this case we have 
non-trivial color structure because the soft gluon makes color connections all the possible 
way between the hard partons. Combining the soft factorization formulae with coUinear one 
our splitting kernel is given by 



,Pm+l)[Vkrn+l{Pk,Pni+l) " Vljn+l{Pl , Pni+l)] 



(6) 



1=1 



k=l 



Here we used the color conservation, that is Tj^ = —Ti ■ J^k^i '^k- One can show that the 
expression under the square brackets vanishes in large angle limit when "dik <C '&im+i, '&km+i 
(The 'dij denotes the angle between momenta pi and Pj.). This effect is know as color 
coherence. Note, the color coherence breaks down with massive hard partons (quarks, SUSY 
particles). In this case we have wide angle soft radiations. 

This is it, we defined a splitting kernel based on the soft and coUinear approximation. 
With this we are able to sum up the leading and next-to-leading logarithms. We still 
have some freedom, for example the definition of splitting kernel away from the limits or 
the momentum mapping. The evolution (ordering) parameter can be basically any infrared 
sensitive variable such as the virtuality or the transverse momentum of the daughter partons. 
Note, if we want to consider spin and color correlations properly in the parton shower we 
cannot avoid negative weights. So far there is only one algorithm has been defined along this 
ideas [S] but it hasn't been implemented yet. From the point of the implementation, the color 
interferences make some complications but one can impose some further approximations to 
simplify it. 

HerwingTD] and old Pythia[TT' implement direct angular ordering [T^. Inserting the 
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conditions for the emission angles into Eq. ([6| then we have 

1=1 k=i 

k^l 

where the dots stand for the neglected terms those are finite in both soft and collinear hmits. 
Note, these contributions are finite only after we perform the integral over the azimuthal 
angle of Pm+i about the direction pi and pk- If we consider only those phase space regions 
where emissions are ordered in angle then 9{'dim+i < i^ik) = 0{-dkm+i < ^ik) = 1 and the 
color part of Eq. ([t]) becomes trivial and the approximated splitting kernel is identical to 
the collinear splitting kernel which is given in Eq. ^. 

The other way to simplify the color structure is to expand the splitting kernel in powers of 
where N,. is the number of the color states in fundamental representation. The gluon 
is a color 8, but in leading color approximation the gluon can be considered to be a 3 (8) 3. 
The partons makes 33 color dipoles with other partons but at leading color level the gluon 
never makes a color dipole with itself. The color connection operator T/ - Tk becomes simple, 
it is non-zero if the partons I and k makes a 33 color dipole and the approximated splitting 
operator is identical to the collinear splitting kernel given in Eq. ([5|. In this approximation 
it is important that the momentum mapping must be exact or based on dipole kinematics; 
if parton I radiates a gluon then the recoiled parton must be the color connected one. 
Ariadne and the new Pythia[11^ implement this approximation. There are some new 
developments |14j based on the leading color approximation and they implement this color 
dipole shower model. 

The parton shower algorithms have been derived from perturbative QCD but we cannot 
consider them as theory predictions because they use rather nonsystematic approximations. 
The original idea was to consider and simulate higher order matrix element by using only 
soft and collinear factorization of the QCD matrix elements. This is a systematical ap- 
proximation since the factorization properties of the matrix elements are held all order. At 
the end of this section it is worthwhile to highlight the addition approximations and the 
limitation of the available parton shower implementations: 

1. The current parton shower programs are still leading order calculations however they 
consider higher order contributions in an approximated way. Dependence on the un- 
physical scales is still strong. 

2. The phase space is usually treated approximately. The angular ordered showers don't 
cover the phase space properly ("dead cone") and some special treatment is required 
to to fill these regions. 

3. The direct angular ordering or the leading color approximation neglect the color cor- 
relations. The color interferences could be significant in the case of non-global observ- 
ables [TS]. Usually the spin correlations are also neglected. Herwig considers spin 
correlations. 

4. They are not defined systematically. The direct angular ordering is not defined or 
hard to define at higher order. Even the kinematics of the dipole shower model is 
inconsistent with the higher order. We have some freedom to define the splitting 
kernel and momentum mapping but the core algorithm should independent of the 
level of the calculation. 
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Figure 5: Calculation of a shower starting with a 2 — > 2 hard cross section (green rounded 
rectangle). The shower evolution operator has been iterated twice, so that the first term 
represents no splitting, the second term has one splitting, the third term has two splittings, 
and the final term contains contributions with three or more splittings. 

5. The only exact matrix element in the calculations is 2 ^ 2 like. If we want to calculate 
say 3, 4, 5, ...-jet cross section we should use 2^3,4, 5, ... LO or NLO matrix elements. 
In the next section I discuss the matching of shower to exact matrix elements. 

6. More questions on non-perturbative effects: What is underlying event? How can we 
model it? How to consider quantum interferences in hadronization models? 

3 Matching parton showers to fix order calculations 
3.1 Born level matching 

The standard shower depicted in Figure |5] has a deficiency, which is illustrated in Figure [6) 
The left-hand picture depicts a term contributing to the standard shower. In this term, 
there are Sudakov factors and 1^2 parton splitting functions. If we omit the Sudakov 
factors, we have the 1^2 parton splittings as depicted in the middle picture. These 
splittings are approximations based on the splitting angles being small or one of the daughter 
partons having small momentum. Thus the shower splitting probability with two splittings 
approximates the exact squared matrix element for 2 — > 4 scattering. The approximation 
is good in parts of the final state phase space, but not in all of it. Thus one might want 
to replace the approximate squared matrix element of the middle picture with the exact 
squared matrix element of the right-hand picture. However, if we use the exact squared 
matrix element, we lack the Sudakov factors. 

One can improve the approximation as illustrated in Figure |7] We reweight the exact 
squared matrix element by the ratio of the shower approximation with Sudakov factors to 
the shower approximation without Sudakov factors. The idea is to insert the Sudakov factors 
into the exact squared matrix element. This is the essential idea in the paper of Catani, 
Krauss, Kuhn, and Webber \W . They use the kx jet algorithm to define the ratio needed 
to calculate the Sudakov reweighting factor. 

There is another way to improve shower as illustrated in Figure [8] First we generate the 
event according to the shower and then rewieght it by the ratio of the exact and approximated 



DIS2007 



Standard shower contribution 



Small pT approximation 











T 


S 








I 



\M{2 ^ 4)|' 



Figure 6: The left-hand picture is the 2^4 cross section in shower approximation. The 
center picture is the shower approximation omitting the Sudakov factors. The right hand 
picture is the exact tree level 2^4 cross section. The cross section based on splitting 
functions (middle picture) is a collinear/soft approximation to this. 





Figure 7: An improved version of the 2 — s- 4 cross section. First we generate the 4-parton 
configuration according to the exact matrix element and take the shower approximation 
(with sudakov factors), divide by the approximate coUinear squared matrix element, and 
multiply by the exact tree level squared matrix element. The graphical symbol on the right 
hand side represents this Sudakov reweighted cross section. 



matrix element. The approximated matrix element is calculated over a unique emission 
history that is determined by a jet algorithm. The original MLM algorithm [17i uses the 
cone algorithm. The advantage of this method over the CKKW method is that the algorithm 
use the native Sudakov factors of the underlying parton shower. 




Figure 8: An improved version of the 2^4 cross section with matrix element reweighting 
factor. 

There is a further step in implementing this idea. CKKW divide the shower evolution 
into two stages, < t < tini and t-mi < t < t{, where ti^i is a parameter that represents 
a moderate Pt scale and t{ represents the very small Pt scale at which showers stop and 
hadronization is simulated. 

With this division, the Sudakov reweighting can be performed for the part of the shower 
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at scale harder than tjni, as depicted in Figure [9] The first term has no sphttings at scale 
harder than ti^i. In the second term there is one splitting, generated via the exact matrix 
element with a Sudakov correction as discussed above. In the next term there are two 
splittings. If we suppose that we do not have exact matrix elements for more than 2^4 
partons, states at scale with more partons are generated with the ordinary parton shower. 
However, this contribution is suppressed by factors of a^. Evolution from iinj to t[ is done 
via the ordinary shower algorithm. 



(i„i tf ti„i tr 




Figure 9: Shower with CKKW jet number matching. The calculation for n jets at scale tiai 
is based on the Sudakov reweighted tree level cross section for the production of n partons. 

To state the main idea of this jet number matching in a little different language, we can 
consider the cross section for an observable F. In the CKKW method, we break (j[F] into 
a sum of contributions ^^[i^] from final states with m jets at resolution scale tini. Then 
crm[-F] is evaluated using the exact tree level matrix element for 2 —^ m parton scattering, 
supplemented by Sudakov reweighting and further supplemented by showering of the m + 2 
partons at scales softer than tini- If F is an infrared safe observable, this method gets 
fm[-F] correct to the leading perturbative order, a™. The method can be extended. The 
present authors have shown (at least for the case of electron-positron annihilation) how to 
get cr,n[F] for an infrared safe observable correct to next-to-leading order, a^~^^ [18]. The 
required NLO adjustments are a little complicated, so I do not discuss them here. 

3.2 Next-to-leading order matching 

Matching parton shower with NLO fix order calculation is a very active field of parton shower 
developments. There are two basic approaches. First one is the MC@NLO project [19]. 
The main idea here is to avoid the double counting by introducing extra counterterm which 
is extracted out from the underlying shower algorithm. This method has been applied for 
several 2 ^ + X where no colored object in the final state and some 2^ l + X, 2 ^ 2 + X 
processes, where the QCD particles in the final state are heavy [TO] . 

The other approach was originally proposed by Kramer and Soper [20J and they imple- 
mented it for e+e^ 3-jets. The idea is to include the first step of the shower in the NLO 
calculation and then start the parton shower from this configuration. Based on this con- 
cept some matching algorithm have been proposed but they are haven't been implemented 
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[131 [IHl |5T] . In the next I discuss in detail only the MC@NLO approach because only this 
scheme has been implemented for LHC processes so far. 

Let us start with the NLO cross section. After applying a subtraction scheme to remove 
the infrared singularities, we have 



aNLo = / [da"" + da"" + da"" + [ da-^jFj") + / 

J rn J 1 J m-\-l 



(8) 



where da^ , da^ , dcr^ , da'^ and da"^ are the Born, real, virtual contributions, coUinear 
counterterm and subtraction term of the NLO scheme, respectively. The physical quantity 
is defined by the functions i^j™^ and iTj™"*"^). 

The naive way to add parton shower corrections is to replace the jet functions with the 
shower interface function. This approach is not good because it leads to double counting. It 
is easy to see, the shower that starts from the Born term generates higher order contributions 
those are already considered by NLO terms. 

To avoid double counting Frixione and Webber P[9J organized the calculation in the 
following way: 



(2^m) 
MC 



/ [da"" + da^ + da^ + f da^]l, 

Jm Jl 

_L f r^^R ^^MC 1 r(2^m+l) / r^^MC J A 1 p 



(2^m) 
MC 



(9) 



Here the contribution da^^^i is extracted from the underlying parton shower algorithm. The 

functions lyl^"^^ and -^mcT™^^' ^'"'^ ^^"^ interface functions to the shower. We have different 
choices for the m and rn + 1 parton interface functions, thus we have 

~Z^(tf,i„) and ~^(if'Wi)A^(Wi,im) • (10) 

In the TO-parton case we simply start the shower from the m-parton configuration while in 
the TO + 1 parton case first we insert some Sudakov factor representing the probability of 
nothing happens between the TO-parton and m + 1 parton states and starts the shower from 
the m + 1 parton configuration. 

There are some limitation of the MC@NLO approach: i) It is worked out for Herwig. 
One has to redo the Monte Carlo subtraction scheme if we want to match say Pythia 
to NLO computations, ii) Matching procedure is defined only for simple processes like 
2 ^ + X, 1 + X, 2. iii) The double counting problem is not fully solved but it is probably 
numerically invisible because of the strong Sudakov suppression. The problem related to the 
soft singularities and it appears only in the 2^2 like processes where the color structure 
is not trivial. 



4 Conclusions 



Parton shower event generators have proved to be an essential tool for particle physics. These 
computer programs perform calculations of cross sections according to an approximation to 
the standard model or some of its possible extensions. Because of the great success of these 
programs, it is worthwhile to investigate possible improvements. 
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In a typical parton shower event generator, the physics is modeled as a process in classical 
statistical mechanics. Some number of partons are produced in a hard interaction. Then 
each parton has a chance to split into two partons, with the probability to split determined 
from an approximation to the theory. Parton splitting continues in this probabilistic style 
until a complete parton shower has developed. 

The underlying approximation is the factorization of amplitudes in the soft or coUinear 
limits. However, further approximations are usually added: i) The interference between a 
diagram in which a soft gluon is emitted from one hard parton and a diagram in which 
the same soft gluon is emitted from another hard parton is treated in an approximate way, 
with the "angular ordering" approximation, ii) Color is treated in an approximate way, 
valid when l/N"^ where A'c = 3 is the number of colors, iii) Parton spin is treated in 
an approximate way. According to the full quantum amplitudes, when a parton splits, the 
angular distribution of the daughter partons depends on the mother parton spin and even on 
the interference between different mother-parton spin states. This dependence is typically 
ignored. With the use of these further approximations, one can get to a formalism in which 
the shower develops according to classical statistical mechanics with a certain evolution 
operator. 

I think the way to improve the parton showers is to formulate it based on the factorization 
of amplitudes in the soft or coUinear limits in which one does not make the additional 
approximations enumerated above. For this, one would have to use quantum statistical 
mechanics instead of classical statistical mechanics. 

On the other hand the parton shower algorithm should cooperate with exact LO and 
NLO matrix elements. Currently we have some very promising tools such as CKKW, MLM 
and MC@NLO matching schemes. The CKKW and MLM matching procedures patch the 
"hole" between the Born level fix order and the shower calculations while the MC@NLO 
and other NLO matching schemes do the same between the shower and fix order NLO 
calculations. If we want more precise tools we need more advanced framework. We need 
a general LO shower framework that naturally includes the LO and NLO calculation. Or 
phrase it differently, we should reformulate the LO and NLO calculation to make the shower 
part of them. 
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